Orbital Lamb shift and mixing of the pseudo-zero-mode Landau levels 

in ABC-stacked trilayer graphene 



(N 

o 

O 

Q 



o3 



O 

o 



> 

oo 
m 
m 
O 

(N 



X 



K. Shizuya 
Yukawa Institute for Theoretical Physics 
Kyoto University, Kyoto 606-8502, Japan 

In a magnetic field graphene trilayers support a characteristic multiplet of 12 zero(-energy)-mode 
Landau levels with a threefold degeneracy in Landau orbitals. It was earlier noted for bilayer 
graphene that Coulombic vacuum fluctuations, specific to graphene, lift the orbital degeneracy of 
such zero-energy modes and that these "Lamb-shfted" orbital modes, with filling, get mixed via the 
Coulomb interaction. It is pointed out that analogous orbital Lamb shift and mixing of zero-mode 
levels can also take place, with an enriched symmetry content, in ABC-stacked trilayer graphene; 
and its consequences are discussed in the light of experimental results. 

PACS numbers: 73.22.Pr,73.43.-f,75.25.Dk 



I. INTRODUCTION 

Graphene, 1-3 an atomic layer of graphite that supports 
massless Dirac fermions, attracts great attention for its 
unique and promising electronic properties. Recently in- 
terest appears to center on bilayers and few layers of 
graphene, where the added layer degree of freedom makes 
the physics and applications of graphene richer. In par- 
ticular, bilayer graphene and some types of multilayers 
enjoy the property that their band gaps are externally 
tunable. 4-7 

A notable signal of Dirac fermions is the fact that 
graphene, in a magnetic field, supports a characteris- 
tic set of zero-energy Landau levels, whose emergence 
and degeneracy have a topological origin in the chiral 
anomaly. 8 Monolayer graphene has four such zero-energy 
levels owing to the spin and valley degeneracy, and they 
are responsible for the observed half-integer quantum 
Hall effect. 1,2 In bilayer graphene there are eight such 
levels, with an extra twofold degeneracy 4 in Landau or- 
bitals n=0 and 1. This "orbital" degeneracy is a con- 
sequence of topology and the added layer, and iV-layer 
graphene necessarily has 47V zero-energy Landau levels 
with TV-fold orbital degeneracy. In the presence of Zee- 
man coupling, Coulomb interactions, etc., these zero- 
energy levels evolve into a variety of pseudo-zero-mode 
(PZM) levels, or broken-symmetry states, as discussed 
theoretically. 9-11 The interplay of orbital degeneracy and 
Coulomb interactions brings about a new realm of quan- 
tum phenomena 9 ' 12-16 in the PZM sector, such as orbital 
mixing and orbital-pseudospin waves. 

Graphene is distinguished from conventional electron 
systems by the feature that it is an intrinsically many- 
body system equipped with the quantum vacuum, or 
the valence band acting as the Dirac sea. Quantum 
fluctuations of the Dirac sea are sizable, even lead- 
ing to ultraviolet divergences; and one encounters such 
field-theoretic (or many-body) phenomena as velocity 
renormalization, 17 screening of charge, 18 and nontriv- 
ial Coulombic corrections to cyclotron resonance. 19-23 
Quantum fluctuations also affect the PZM levels sub- 



stantially. They work to lift 24 the orbital degeneracy of 
the PZM levels in bilayer graphene; each orbital mode 
responds to quantum fluctuations differently and gets 
shifted, just like the Lamb shift 25 in the hydrogen atom, 
where the field-theoretic effect of quantum electrodynam- 
ics was revealed for the first time historically. The Lamb- 
shifted orbital modes get mixed via Coulomb interactions 
and govern the fine structure of the PZM sector. 

A number of recent experiments 26-32 have verified that 
the electronic properties of graphene trilayers strongly 
depend on the stacking order, with Bernal (ABA)- 
stacked trilayers remaining metallic in contrast to rhom- 
bohedral (ABC-stacked) trilayers which exhibit a tun- 
able band gap. Actually trilayers drew theorists' at- 
tention 33-36 even before experiments and their rich elec- 
tronic properties 37-39 have been under active study. Cur- 
rently considerable attention 40-42 is directed to ABC- 
stacked trilayers which are a chiral generalization of bi- 
layer graphene. In view of this, it is of interest to ask how 
the Coulombic vacuum and orbital dynamics generalizes 
to trilayers. 

The purpose of this paper is to examine the effect 
of Coulombic vacuum fluctuations in trilayers and show 
that the orbital Lamb shift and orbital mixing of the 
PZM levels are also present, with an enriched symmetry 
content, in A73C-stacked trilayer graphene. It is noted, 
in particular, that level mixing takes place without level 
crossing; this mechanism would, for high-quality sam- 
ples, lead to an observable sequence of fully-split broken- 
symmetry quantum Hall states in the PZM sector. 

In Sec. II we briefly review some basic features of the 
PZM levels in ABC-trilayer graphene, and in Sec. Ill 
show that vacuum fluctuations lift their orbital degener- 
acy. In Sec. IV we discuss in a simplified setting how 
orbital mixing of the PZM levels takes place via the 
Coulomb interaction. In Sec. V we examine the hier- 
archy of broken-symmetry states under practical condi- 
tions. Section VI is devoted to a summary and discussion. 
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II. TRILAYER GRAPHENE 

The A_BC-stacked trilayer graphene consists of three 
graphene layers with vertically-arranged dimer bonds 
(Bi,A 2 ) and (B 2 ,A 3 ), where (A^Bi) denote inequiva- 
lent lattice sites in the i-th layer. The interlayer cou- 
pling 7q = jBiAi ~ 3 eV is related to the Fermi velocity 
v = (\/3/2)aL7o/?l ~ 10 6 m/s in monolayer graphene. 
Interlayer hopping via the nearest-neighbor dimer cou- 
pling 43 71 = -fB 1 A 2 — 1b 2 a 3 ~ 0.4 eV leads to soft cubic 
spectra 33 oc |p| 3 in the low-energy branches |e| < 71. 

The effective Hamiltonian for ABC-stacked trilayer 
graphene with such intralayer and interlayer couplings 
is written as 36 



iP" = 

Uk = 

A = 
W = 




Hk* k + {^ k ') ] U K '^ 



K 



V = 



v 4 pt V3 p 



7i 



— V4 pi 



(1) 



with p — p x + ip y and pi = p x — ip y . Here ^> K — 
(ij)A 1 ,ipB 1 ,ipAi,i>B2i' l l ) A 3 ,ipB 3 ) t stands for the electron 
field at the K valley. V3 and W4 are related to the non- 
leading interlayer couplings 73 = r jA 1 B 2 and 74 = Ja 1 a 2 , 
respectively, and 72 = 7AiB 3 - (^1,^2,^3) stand for the 
on-site energies of the three layers; we take U 2 — with- 
out loss of generality. As in bilayer graphene, 4 these bi- 
ases {Ui} open a tunable band gap 33 ~ U\ — U3. Hk is 
diagonal in (suppressed) electron spin. 

The Hamiltonian Hk' at another valley is given by Hk 



with p 



-Px + ip v 



—p^ and pi — > — p, and acts on a 



spinor of the same sublattice content as f x . Actually, 
Hk' is unitarily equivalent to Hk with the sign of V3 and 
72 reversed and with layer 1 and layer 3 interchanged, 



S*Uk>S 



h 



K |-u 3 ,-72;£/i<->f/3> 



02 



-cr 2 



(2) 
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In view of this, we adopt Hk 1 — S^T-Lk'S for Hk 1 
and simply pass to the K' valley by reversing the sign 
of V3 and 72 and interchanging U\ and U3 in the K- 
valley expressions. Nonzero v 3 , 72 and bias U\ — U3 
thus act as valley-symmetry breakings. Remember that 
in this representation Hk 1 acts on a spinor of the form 

= (lpB 3 , VU 3 , *Pb 2 ,1pA 2 ,1pB 1 ,1pA 1 ) t - 

A direct link between the K- and if'-valley represen- 
tations, such as Eq. (2), was also noted 24 for bilayer 
graphene. We remark that such a link is not shared by 
ABA-stacked trilayers, where the Landau-level spectra 
significantly differ 37 between the two valleys for nonzero 
biases (though they coincide for zero bias). 



We have discussed the general structure of trilayer pa- 
rameters for completeness. For our present analysis of 
quantum effects in ABC trilayers we retain only the lead- 
ing parameters (v,-fi,Ui); the effect of nonleading cou- 
plings («3, V4, 72) is discussed later in Sec. V. In addition, 
we focus on the case of a symmetric bias 33 by choosing 
U 3 = -U x = u/2. 

Let us place trilayer graphene in a strong uniform mag- 
netic field B z — B > normal to the sample plane; we 
set, in Hk, P — > n = p + eA with A = A x + iA y = —B y, 



and denote the the magnetic length as i = 1/veB; set- 
ting a = \lleB IT+ then yields [a, a T ] = 1. It is easily seen 
that the eigenmodes of Hk have the structure 



(|n-3) bW,\n-2)dg\\n-2) b™, 
In-l)^),^-!)^),^)^))* 



(3) 



with n = 0, 1, 2, ■ ■ ■, where only the orbital eigenmodes 
are shown using the standard harmonic-oscillator basis 
{\n)} (with the understanding that \n) = for n < 0). 

The Coefficients V n = (b i n\d { n ) ,b { n ) ,d i n\b i n\d i n ) ) t for 

each n = 3, 4, . . . are given by the eigenvectors (chosen to 
form an orthonormal basis) of the reduced Hamiltonian 
"Hrcd = u c H n with 



/ —M 

y/n-2 



H n — 



where 



V^2 
-M 

7 



\ 



7 






7 



7 
M 



s/n M ) 



(4) 



to c = V2v/£ w 36.3 x u[10 6 m/s] y/B[T] meV, (5) 

with v in units of 10 6 m/s and B in tesla, is the charac- 
teristic cyclotron energy for monolayer graphene; M = 
2 u/lj c and 7 = 7i/w c . Note that eigenvectors v„ can be 
taken real since H n is a real symmetric matrix. 

Solving the secular equation shows that there are 6 
branches of Landau levels for each integer n > 3. We 
denote the eigenvalues as e_„» < e_„/ < e_„ < < e„ < 
£«' < so that the index ±n reflects the sign of e„; 
|e±n'| ^ 7i and |e± n «| > 71. The \n\ = 3 levels, e.g., 
consist of the n = (±3. ± 3', ±3") branches. 

There are also solutions for n =2, 1 and 0, for which 
H n is reduced to a matrix of smaller rank 5, 3 and 1. For 
n = 0, Hred has an obvious eigenvalue eo = t/3 = u/2 
with eigenvector v = (0, 0, 0, 0, 0, 1)* or 



*o = (0,0, 0,0, 0,10))*. 



(6) 



For n=l, "Hred has three eigenvalues (ei,e±i')> which, 
for u = 0, read (0, ±-\/7 2 + l)ui c . The zero-energy solu- 
tion, in particular, takes the form: 



*i Cl (0,0,0, 



-^10>,0, |1>)*, 



(7) 
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with k = I/7 and c\ — j/\/j 2 + 1 = + k 2 . For 

n = 2, 7i rcc i has five eigenvalues (e 2 , e± 2 <, £±2"), with 
e 2 = and |e± 2 '| ~ |e±2"| ~ 7i for u = 0. The zero- 
energy solution takes the form 

* 2 u =° c 2 (0,V2k 2 |0),0, - V2 k 0, 12))', (8) 

with c 2 = 7 2 /^2 + 2f + 7 4 = 1/V1 + 2k 2 + 2k 4 . Note 
that these zero-energy solutions (^Fo, ^2) reside pre- 
dominantly on the i? 3 lattice sites of the third layer; cor- 
respondingly, the zero-energy solutions at the if' valley 
reside predominantly on the A\ sites of the first layer. 

Of our particular concern are these three zero-energy 
modes (^o, ^2)- For u = there are 12 such zero- 
energy Landau levels differing in spin, valley and orbital 
[n = (0,1,2)] degrees of freedom; their presence is dic- 
tated by the nonzero index 8 ' 18 of the Dirac Hamiltonian 
Hk © W.K' with only v and 71 retained. 

For nonzero bias 11/O they evolve into the pseudo- 
zero modes with nonzero energies, 

(e%,e^,e 2 i ) = (l,l-z 1 ,l-z 2 )u/2, 
Zl = n 2 { Cl ) 2 +0(u 2 k 6 ), 

z 2 = 2n 2 {l + 2n 2 ){c 2 ) 2 +0{u 2 n i ) 1 (9) 

where u = u/uj c . One can also write e" ~ (ci) 2 u/2 and 
e^^ (1-2k 4 )(c 2 ) 2 u/2. 

For a numerical estimate let us take, as typical val- 
ues, 43 70 = 3.16 eV (or v w 1.0 x 10 6 m/s) and 71 = 0.4 
cV. They yield 7 w 3.41 and k w 0.293 at B = 10T, 
which in turn lead to c x w 0.960, c 2 w 0.918, z x w 0.08 
and z 2 w 0.17. One thus has 

(eg,ey,e^)w (1,0.92, 0.83) u/2 (10) 

for u -C co c . 

One can pass to the if' valley by setting u — » — u in 
the if- valley expressions. The eigensy stems (e„,v„) at 
the two valleys are related as 

£n\K' = — £-n\K, 

bW\K' = -b { \\ K ,d$\ K ' = d^ n \ K (11) 

for each mode (n,n',n") and i e (1,2,3). For later con- 
venience, we continue to use n — (0, 1,2) to specify the 
PZM levels at the if' valley; one can thus effectively set 
n = ±0 -> 0, ±1 -> 1 and ±2 -> 2. When the interlayer 
bias u is turned on, these PZM levels go up or down 
oppositely at the two valleys, opening a band gap <~ u. 

The Landau-level structure is made explicit by pass- 
ing to the \n, yo) basis (with yo = i 2 p x ) via the ex- 
pansion (**(x),**'(x)) = E„,„ yo) {C ;Q (l/o)}, 
where n refers to the level index, a e (j, I) to the 
spin, and a <G (if, if') to the valley. The charge density 
P-p = Jd 2 xe l P x p with p = (^ K )^ K + (^')t^' i s 
thereby written as 24 

00 

\ " \ " kn;a r>kn;aa 

H—p — IP / j / j 9p ll aa;p ) 

fc,n= — 00 a, a 



<Tp = Jdyo€' af (yo)e^^ b (y ), (12) 

where 7 P = e~ f2p2 / 4 ; r = (i£ 2 d/dyo,yo) stands for the 
center coordinate with uncertainty [r x ,r y ] = il 2 ; the level 
sum is taken over possible (n, n', n"). 

The coefficient matrix g^ a = 3p"| a at valley a £ 
(if, if') is constructed from the eigenvectors v n | OJ 

s*» = 6«6£>/i*i-».n-» 

+(4 1) 4 1) +4 2) C ) )/if | - 2 ' M - 2 
+(4 2) 4 2) +4 3) ^ 3) )/if | - 1 ' M - 1 
+4 3) 4 3) / P fc|>l , (13) 

where 

/ p fc " = v^!7fc!(- g 7y2) fc -"4 fc -")(|g| 2 /2) (14) 

for k > n > 0, and /£* = (/ fo p )t ; g = £( Px - lPy ) ; it is 
understood that = for fc < or n < 0. As seen 
from Eq. (11), <7p" ;a at the two valleys are related as 

rt mn| -,— m,— roi „mn:a\ „—m,—n\a\ {-\K\ 

9 P \k> — 9 P \k, 9 P \u — 9 P \-u- (toj 

Within the n e (0, 1, 2) sector, g£" ;a are functions 44 of 
(u 2 , k 2 ) and are thus common to both valleys; for u = 0, 
they read 

3p = 1, = Cl £p/V2, g p ° = - Cl £p/V2, 
g° p 2 = c 2 £V/(2V2), 
gp 1 = 1-( C1 ) 2 KP 2 , . 9p 2 = Ap. 9 ° 4 , 
4 2 = l-(c 2 ) 2 [(l + K 2 )^ 2 p 2 -i(£ 2 p 2 ) 2 ], 
A p = V2(l + « 2 -i£ 2 p 2 ) C2 , (16) 

with ci = 1/VI + k 2 and c 2 = 1/V1 + 2k 2 + 2k 4 . 

From now on we frequently suppress summations over 
levels n, spins a and valleys a, with the convention that 
the sum is taken over repeated indices. The Hamiltonian 
H hl projected to the PZM levels is thereby written as 

H u = e u n SRpp.0 - fiz {T-i)0 a R^sfo (17) 

with n e (0,1,2) and SR^ = R^ K - R^f*'. 
Here the Zeeman term = g* urB « 0.12 B[T] meV is 
introduced via the spin matrix T3 =0-3/2. 

III. VACUUM FLUCTUATIONS 

In this section we examine the effect of Coulom- 
bic quantum fluctuations on the PZM multiplet. The 
Coulomb interaction is written as 

V = \^2 v p -P-pPp : > (18) 
p 
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where Dp = 2ira/ (eb|p|) with a = e 2 / (47reo) ~ 1/137 and 
the substrate dielectric constant e^; ^ p = J cPp/(27r) 2 . 
For simplicity we ignore the difference between the in- 
tralayer and interlayer Coulomb potentials. 

In this paper we generally focus on many-body ground 
states \G) with a homogeneous density, realized at integer 
filling factor v G [—6,6]. We set the expectation values 
(G\R™g b \G) = S KoP0 ^p ab with p Q = l/(2^ 2 ) and 
^k,o = (27r) 2 <5 2 (k); accordingly, the filling factor v^ aa = 
1 for a filled level specified by (n, a, a). 

Let us define the Dirac sea |DS) as the valence band 
with levels below the PZM sector (i.e., levels with n < 
-3, n' < -1' and n" < -2") all filled. We construct the 
Hartree-Fock Hamiltonian V HF out of V as the effective 
Hamiltonian that governs the electron states over |DS). 
Let us write V HF = Vb + Vx ■ As usual, the direct inter- 
action Vd oc Vp^ R™J!^ ' is removed if one takes into 
account neutralizing positive background charges. We 
thus focus on the exchange interaction 

\r — \^ ,, .,2 mn';b m'n;a , mn;ba jjm'n'iab / 10 \ 

p 

where we sum over filled levels (to, n) and retain the PZM 
sector to', n' G (0, 1, 2). 

Let us first extract, out of Vx, the contribution from 
the Dirac sea, 

t/DS _ \ .. _,2 \ " I m'n;a\2 n™' W ' ;aa / on \ 
V X --2^ V P^P 1^ l5p I H <*a;0 > ( 2 0) 
p riGDS 

where the sum over to' G (0, 1,2), a G (K,K') and a G 
(t, |) is understood. Actually, the sum over infinitely 
many filled levels with — oo < n G DS gives rise to an 
ultraviolet divergence. 

Fortunately one can isolate the divergence and even 
evaluate V^ S exactly for zero bias u — »■ 0, as done for 
the bilayer case. 24 Note first that, as seen from Eq. (15), 
9p n \K = 9p n \K' = Sp m '~™ for u = 0, and use the com- 
pleteness relation 24 

oo 

E l5p"| 2 = e£V/2 ( 21 ) 

n=— oo 

to extend the sum X)«eDS to complement X^eDS as 
well. The result is 

E K\ 2 U =\^ 2/2 -\9ir-\9^ 2 -W v 2 \% (22) 
neDS 

for j G (0, 1,2). Equation (21) was noted earlier with a 
formal proof; a direct proof of it is given in Appendix A. 
The e e p / 2 term in Eq. (22), though leading to a diver- 
gence upon integration over p, is common to all levels 
j and is safely omitted. We thus take the rest as the 
regularized expression for XneDS lsp™l 2 - 

The regularized Dirac-sea contribution thus reads 

l/DS v D 00;aa . v u ll;aa . v u 22;aa (nn\ 

V X = e H aa;0 + e l H aa;0 + e 2 H aa ;0 i 12JJ 



«J = *Ev£El4 B | a - (24) 

p n=0 

Integration over p, with the aid of the formula 

E«p7p[l,<zW,<zV] = [1,1, 3, 15, 105] V c (25) 
p 

with q = £\p\, then yields 

4 = Hi + c2c i(l + f- 2 + « 4 )]^, 

el = i [1 + cfc 2 (jj + if k 2 + k*)]V c , 

e 2 = I [l + c l c 2 (|| — 64 r2 — "F ft4 _ _ 2k 8 )] V c , 

(26) 

where c 2 = (ci) 2 , etc., V c = y^n~/2V c and 

K = a/(e b £) « (56.1/e b ) /B[T] meV. (27) 
Numerically, 

(eg, e^, e£) = (0.888, 0.777, 0.641) V c (28) 

for 7 = 1/k w 3.41 at S = 10 T. 

Vacuum fluctuations thus shift the n = 0,1 and 2 
modes differently and the splitting among (eoi e i\ e 2) re ~ 
fleets the difference in their spatial distributions, as is 
clear from Eq. (24). The empty PZM levels are ordered 
as 6q > e\ > e\ > 0. Actually the spectra vary with 
filling of the PZM sector. Note Eq. (20), which tells us 
to include extra contributions — |g-'"| 2 for ej, when the 
n G (0, 1, 2) level is filled. In particular, when the PZM 
sector is filled up, one finds that {ej} change sign so that 
they are ordered as Cq < e\ < e\ < 0. 

Let us next suppose filling the lowest-lying n = 2 level 
first in the empty PZM sector (for u = p z =0). One 
then finds t\ w 0.054 V c for g w 3.41. If, instead, the 
highest-lying n — level were first filled, one would find 
eg rj —0.11 V c . This puzzling situation suggests that one 
cannot reach the true ground state by filling the n = 2 
level alone. It is clear now that one has to diagonalize 
the exchange interaction (19), with mixing among the 
n = (0, 1, 2) orbital modes taken into account. 

IV. MIXING OF THE PZM LEVELS 

In this section we examine how the PZM sector changes 
in spectrum with filling. The first step is to extract from 
Vx in Eq. (19) the exchange interaction acting within the 
n = (0, 1, 2) sector, 

x = -Z^P^p r p ' 

P 

pOO = ^nO^rll = ^^l^pffl = ^«| ff «2 |2j 

C = 5pV> 10 + . 9p 1 3 2 -> 21 , 
rp 1 - ^p^ + ffp 1 ^^ 21 , 
r 20 = .9p°g 22 P ^ 20 , (29) 
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where m,n £ (0,1,2); T" 1 = (r 1 ^)*, r 21 = (rL 2 p )*, etc. 
For conciseness the spin and valley indices have been sup- 
pressed in the above; v nn R °, e.g., reads ^™' &a R^p^o- 

Let us, for the moment, freeze the spin and valley 
degrees of freedom and focus on the orbital degrees of 
freedom. The PZM sector then consists of three lev- 
els n £ (0,1,2) governed by the effective Hamiltonian 
V = V£ s + Vf = H mn R^ n with H mn = e v n 5 rnn - 

E P u p7p r p™- Note that r p " are real for real nllin g Ac- 
tors v mn , which we take. It therefore suffices to use a 
real 0(3) rotation, rather than a full SU(3) rotation, to 
diagonalize the 3x3 real symmetric matrix H mn . We 
thus rotate ip m = (ip°, ip , -ip 2 ) in orbital space, 

^ m (yo) = [U(e 2 ,9 1 ,e )} mn ^ n (y ), (30) 

with three Euler angles (9 2 ,9±,9o) parameterizing 

U(9 2 ,9 1 ,9 ) = e ie 2 * V Sl * V e ° *°, (31) 

where the spin-1 generators (t a ) bc = ie bac in terms of the 
totally antisymmetric tensor e abc with e 012 = 1. Note 
that 6> mixes n = (1,2), 9\ mixes (0, 2), etc. 

Via the rotation, V = H mn Bg n = U mn TZ^ n with H = 
WHU, where lZ nn stand for the charge operators for 
i.e., R™p.Q with ipp' — > . The transformed 
fields <E>™ are taken to diagonalize H mn and hence the 
associated filling factors as well, N n oc (G\(^ n )^ n \G) 
with < N n < 1 and n £ (0, 1,2); one can now write 

V mn = (CW)*JVn'(^*)^n. 

Let us start filling the empty PZM sector at (relative) 
filling factor rif — 0. Obviously, in view of level split- 
ting (28), it is the lowest-lying n — 2 level ($ 2 ) that 
starts to be filled. To follow how it evolves let us sup- 
pose that it is filled with fraction nf < 1 and substitute 
(N , Ni,N 2 ) = (0, 0, n f ). H mn is diagonalized if one can 
adjust (6»o, 0i, 6 2 ) so that H 01 = H 02 = H 12 = 0. 

Note first that, with no level mixing, i.e., do = 9% = 
62 = 0, the eigenvalues {H nn } simply go down with in- 
creasing nf. Note next that, to first order in {0„}, 

H 12 « (0.136- 0.268 n f ) O + ..., 

n 02 - (0.247- 0.130 ra f ) 0i + ..., 

U m ks -0.201 n f 6»o + (0.111 + 0.0626 n f ) 2 + ... (32) 

This structure reveals that 0o = 01 = 02 = for ri{ < 
n cr ~ 0.507 while 9q =/= is possible for nf > n cr . Solving 
for {9 n } numerically for nf > n cr shows that the energy 
eigenvalue H 22 is indeed lowered for nf > n cr with 9q ^ 
0. One can then reach the nf = 1 state, and setting 
(No, Ni, N 2 ) — > (0, nf — 1,1), etc., takes one further to 
the nf =2 and 3 states. 

Figure 1 (a) shows how angles {9 n } vary as nf is in- 
creased from to 3. Actually we find another solu- 
tion which differs from one shown in the figure by signs, 
(0 2 , 0i, 0o ) — > (~02, 0i, — 0o)- These two solutions are re- 
lated by a unitary transformation Y — diag[— 1, 1, — 1], 
withW(-0 2 ,0i,-0 o ) = YU(9 2l 9 1 ,9 )Y~ 1 , so that Yip = 





FIG. 1: Orbital mixing, (a) Angles (60,6-1,62) vary from 
to 7r/2 with filling of the ($°,$ 1 ,$ 2 ) sector, (b) Variations 
of the spectra of the and $ 2 levels over the range of 

relative filling factor nf £ [0,3]. (c) Electric dipole moment 
induced via orbital mixing. 



U(— 02, 0i, — 0o)V$. They naturally lead to the same 
level spectra {e$r,} depicted in Fig. 1 (b). 

In Fig. 1 (a), each 0„ evolves from to 7r/2 with in- 
creasing nf. The eigenmodes ($°, $ 1 , $ 2 ) thereby contin- 
uously change from (tp°, ip , ip 2 ) to (ip 2 , — ip 1 , ip°). The 
empty n = (0, 1,2) levels at rif = thus turn into the 
filled n = (2,1,0) levels of energies (-|e£|, -|e^|, -|e^|), 
respectively, at nf = 3, without any level crossing. Each 
spectrum e$n (nf ) goes down with nf , with marked change 
across nf ~ (0.51, 1, 2, 2.49). The spectra as a whole real- 
ize particle-hole symmetry, with e$o(nf) = — e$2(3 — nf) 
and e$i(nf) = — e$i(3 — nf) in obvious notation. In par- 
ticular, the band gaps at nf =1 and 2 are equal, with 

e gap U=i,2«0.53V c , (33) 

considerably smaller than the full Coulombic gap 2 e 2 « 
1.28 V;. 

A special feature associated with orbital mixing is that 
charge carriers acquire electric dipole moment, as noted 
earlier 12,14 for bilayer graphene. To see this let us con- 
sider coupling to an external scalar potential Aq, with 
the Hamiltonian Ha — ~eE P (^o)p+p- Note that 
Sp 1 , <7p 2 (in P-p) oc p, which implies that orbital mix- 
ing gives rise to coupling to an inplane electric field 
E|| = (E x ,E y ) = —\7A . Indeed, for a spatially almost 
uniform field E|| , the relevant portion of Ha is written as 
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H A a h% n R™ with h A = -oa (ei/V5) (S a T y + £ X T X ), 
where 



(i) < M < M c 



(ii) U cr < U 




act on fields (tp° , ip 1 , i/; 2 ) 1 and A ee A p=0 = V2(1 + k 2 )c 2 . 
The expectation value (G\H A \G) then reads 



(G\H A \G)^ Po yd 2 x(-d-E||), 



(35) 



where d = (0,d v ) and d y — (e£/y/2)ciD(ri{) with func- 
tion D(n { ) = {T y + \T x ) mn v mn given by the plot in Fig 1 
(c) . This shows that electrons acquire electric dipole mo- 
ment of magnitude |d e | = (e£/v2) ci\D(n pz )\ (per parti- 
cle) , pointing in the y axis for the present choice of wave 
functions. 

Actually the electric dipole can point in any direction 
(in general, perpendicular to traveling waves) at no cost 
of energy. To see this let us consider a phase rotation 
of the form (within SU(3) rotations), ip n — > ip n with 
(iP ,^ 1 ,^ 2 ) = (e- i ^-4>°,ip 1 ,e i ' t '-4> 2 ). Note that g* n are 
thereby transformed so that /3_ p remains invariant. This 
transformation leaves V = V® + V£ z invariant; hence, 
the spectrum remains unchanged. Still the electric dipole 
thereby rotates so that 



d = (d x ,d y ) = (sin 0, cos 0) |d e 



(36) 



It is now clear that a pair of solutions tp and Yip, encoun- 
tered earlier, differ by a rotation by tt of coordinates in 
the sample plane. 



V. GENERALIZATION 

In this section we recover the electron spin and valley 
degrees of freedom and explore the PZM sector with both 
Hz and bias u, using the full Hamiltonian 



V 



DS 



v x ■ 



(37) 



We leave u arbitrary but keep \u\ <C V c so that one can 
still use the u = expressions for V = Vx S + V^ z , with 
u retained only in H u as a small perturbation. 45 

In addition, we ignore the difference between the intra- 
and interlayer Coulomb potentials that leads to a valley- 
symmetry breaking of 0(V c d/£), with the layer separa- 
tion d ~ 0.35 nm -C I. This breaking contains capaci- 
tance energies that determine how valleys rotate. In con- 
ventional bilayer systems, the capacitance energy, though 
as tiny as 0(V C d 2 /£ 2 ), is positive and induces a valley ro- 
tation (K,K') — > K ± K 1 , which makes the symmetric 
states oc K + K' lower in energy. In contrast, for bilayer 
graphene, capacitance energies turn out to be negative 24 
and suppress possible valley rotations for u ~ 0. 

Experimentally, it is difficult to directly observe valley 
quantum numbers, especially from the sequence in which 



■t) =f -_^± 



i-i _L±i 



el — - 



0. 



1-T 1+ ' 



1 +T i 



l-i 



2-- 



FIG. 2: Empty levels in the PZM sector at v = -6; for 
illustration, (i) u — 0.25 u CI and (ii) u = lAu cr with u C r = 
/Uz/(1 - z 2 ) ~ 1.2 ^ z and (iz/Vc = 0.05. 



the broken-symmetry states emerge with varying filling 
factor v or magnetic field B. The sequence is governed 
by the Coulombic gaps, which, though possibly triggered 
by small valley or spin or orbital breaking, are practically 
insensitive in magnitude to small |u| <C V c . (In contrast, 
for large bias u, the valley is naturally polarized in either 
K or K' , depending on the sign of u.) For this reason, 
instead of a (rather laborious) analysis of capacitance en- 
ergies, we here simply suppose a possible valley rotation 
(K, K 1 ) — > (+, — ) without specifying its details for small 
u; we take the (-) state to be lower in energy for each 
n G (0,1,2) and spin a G (t,4-)- 

In H e e the exchange interaction V® s + V^ z conserves 
both valley and spin, but breaks the orbital degeneracy. 
In contrast, the small perturbation H u lifts all three de- 
generacies. Figure 2 depicts the empty PZM sector (at 
v = —6) governed by H u + V^ s , with level spectra 

e^ = e v n ±el-y z , e^ = e v n ±e^ + ^ z , (38) 

in obvious notation. There are two possible level pat- 
terns, depending on (i) < u < u CT (of spin-breaking 
domination) or (ii) u > u CT (of valley-breaking domina- 
tion) withu cr = ^ z /(l-z 2 ) [ra 1.2 Hz at 10T]. In Eq. (38), 
for definiteness, we have assumed no valley rotation and 
u > 0, so that (+, — ) = (K,K r ). When a valley rotation 
is induced (for small u), the ±e" portions are replaced by 
more complex expressions, which, at any rate, are small 
for small u, and the level pattern (i) in Fig. 2 remains 
essentially intact. (For consistency, we set u —¥ in our 
discussion for case (i) below.) 

Let us start filling the empty PZM sector. Obviously, 
with \u\,/j,z <C V c , the orbital splitting among {e^} sin- 
gles out the 2_-j- level as the lowest-lying one in both cases 
(i) and (ii). It is thus the 2_-j- level that is filled first. As 
it is being filled, it comes down in energy, followed by 
the 0_^ and l_f levels coupled via the exchange inter- 
action Vx '. These three levels undergo orbital mixing, 
discussed in the previous section, through the v = — 5 
and -4 states until one reaches the v = — 3 state, which 
is orbitally neutral (an SU (3) singlet) but is polarized in 
valley and spin (— , f). The associated v = — 3 level gap 
is a valley gap for case (i) and a spin gap for case (ii), 



£ iy = -3l(l) 

e i/=-3l(n) 



2^, 



(39) 
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FIG. 3: Spectra of the PZM Landau levels at each integer filling factor v £ [—6,6]. (a) u = 0.4 /iz < u C i and (b) u — 
2 nz > u CI with Hz/Vc = 0.05 and u cr = /iz/(l — £2) ~ 1-2 /iz for illustration. Large blobs and squares refer to spin-up levels 

(2_-|-, -|- ) I e and (2+f , 1+f, 0+-f-)|e, respectively, from left to right; smaller symbols refer to those for the spin-down levels. 

Here (2_-j-, 0_-j-)|g, e.g., stands for (2_-|-, 0_-f- ) for empty levels (8 = 0) and (0_-|-, 2_^) for filled levels (8 — 7r/2), 
in accordance with Fig. 1; empty levels have positive energy and occupied levels have negative energy. Valley indices ± arc 
attached to some symbols to indicate the nature of the associated gaps. 



with 2 €\ ~ 1.28 V c . Similarly, as one goes up from 
v = — 3 to = 0, essentially the same orbital mixing 
is repeated for the (+,t) sector in case (i) and for the 
(— ,1) sector in case (ii); analogously for the v 6 [0,6] 
domain. 

Figures 3 (a) and 3 (b) show the resulting spectra of the 
PZM multiplet at each integer filling factor v g [—6,6]. 
They differ in pattern for (i) u ~ and (ii) u > u CI , 
but form a perfectly particle-hole symmetric spectrum 
for the PZM sector in each case. The v — ±2,±1 and 
states thus differ in composition, depending on u. The 
v = state, in particular, is spin-polarized for u ~ and 
valley-polarized for u > u CI , with a gap 

4 a = P ol(i) « + 

4=o I (ii) = 2el+(l-z 2 )u-fi z . (40) 

As to the v = — 5 gap, especially for case (ii) (of rela- 
tively large u), we note the following: (€^,€1,62) in H u , 
via the rotation U, turns into (0.96, 0.899, 0.897) u/2, i.e., 
e" sa €2 at v = — 5; similarly, £q « e % { at v = —4. This 
suggests that the v = (±5, ±4, ±2, ±1) gaps are prac- 
tically insensitive to both bias u and /xz, and equal to 
e gap k=i,2 in Eq. (33), 

4 a = P ±5,±4, ± 2,±i« 0.53 V;. (41) 

These orbital gaps are considerably smaller than the 
(Coulomb-enhanced v = 0, ±3) spin or valley gaps, 

,gap < gap < gap / gap x / 

in conformity with Hund's rule. 9,40 These v = ±1, ±2, ... 
orbital gaps and the v = ±3 valley gaps for u ~ barely 



depend on fiz and will therefore be insensitive to an addi- 
tional parallel field Bu in experiments with a tilted mag- 
netic field, in contrast to the v — spin gap for u ~ 0. 
For bilayer graphene, the corresponding gaps take place 
at v = — 3 and -2, and it was observed 46 that the associ- 
ated resistance minima are barely affected by Bu . 

The orbitally polarized states at v — ±1,±2,±4, ±5 
have spontaneous electric dipole moment and may poten- 
tially be unstable 12,41 ' 42 against charge inhomogeneities. 
Their spectra may be modified (in random patterns or 
regular 41,42 patterns) around local charge concentrations 
but, as long as the orbital gaps survive, the quantum 
Hall states would emerge. Such an instability disappears 
when bias u is sufficiently large to stabilize the valley- 
polarized states. For bilayer graphene full splitting of 
the PZM levels has indeed been observed. 46,47 

The transport properties of trilayers have been studied 
in a number of experiments. 26-32 Experimentally there is 
clear evidence for formation of the quantum Hall states 
in the basic filling-factor sequence v — ±4(iV + 3/2) = 
±6, ±10, ±14, ... for both ABC-and ABA-stacked trilay- 
ers. Evidence is yet very limited for the fine structure 
of the PZM sector with \v\ < 6 in ABC trilayers: An 
experiment, 31 using a Hall-bar device, observed a weak 
anomaly in a xy indicative of the developing v = ±3 gap. 
A clear signal for the v — gap comes from the observa- 
tion 26,30 of the insulating state at the Dirac point (v = 0) 
in ABC-trilayer devices, both suspended and substrate- 
supported ones, with the resistance rising exponentially 
with increasing B and lowering temperature T. Experi- 
mentally, it is normally the v = insulating state that is 
first observed as a nontrivial feature within the PZM sec- 
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tor of few-layer graphene. This suggests that the v = 
gap is an interaction-enhanced gap rather than the far 
smaller intrinsic spin or valley gap. The v = ±3 gaps 
will be the next to be visible via quantized conductance. 
In view of Eqs. (39) and (40), the v = gap will become 
even more prominent with increasing bias u, in contrast 
to the v — ±3 gap. 

Finally we wish to discuss possible effects of nonlead- 
ing interlayer couplings (1*4,^3, 72)- The effect of v 4 can 
be included in % n of Eq. (4) while V3 and 72 induce 
transitions that go outside the PZM sector, as seen from 
the solutions in Eqs. (6) - (8). Accordingly the spectra 
( e 0i e ii £2) are corrected to first order in v 4 /v and to sec- 
ond order in v 3 /v and 72/70- With typical values 28 ' 43 
v 4 /v = r 4 ~ 0.01, vi/v — 73/70 ~ 0.1, and 72 <~ —0.02 
eV, such corrections are generally small. The leading 
0(114) corrections, in particular, may conveniently be in- 
cluded in H u if one sets €q — u/2, e" = (1 — z\)(u/2 + 
2nr 4 uo c ) and eJf = (1 — 2 2 )m/2 + 4/t(l + n 2 )r 4 uj c , with 
2kt- 4 ~ 0.005 and 4k(1 + n 2 )r 4 ~ 0.01. Unlike u, such 
0{vi) corrections are common to the K and K' valleys 
and lead to weak electron-hole asymmetry. The rela- 
tive magnitude of (eg\ e i! e 2) mav vary with bias u and 
can potentially control a valley rotation for small u ~ 0. 
Still the orbital splitting among {e^} is generally larger 
than the splitting among {e^}, and the PZM sector will 
essentially maintain the spectra shown in Fig. 3. The 
electron-hole symmetric spectra there will also serve as 
the base point for further examining possible effects of 
nonleading intra- and interlayer parameters. 

VI. SUMMARY AND DISCUSSION 

In a magnetic field graphene trilayers acquire, on topo- 
logical grounds, a special multiplet of nearly-zero-energy 
Landau levels with a threefold degeneracy in Landau or- 
bitals. In this paper we have studied the structure of 
this PZM multiplet in ABC-stacked trilayer graphene 
and pointed out that its orbital degeneracy is lifted by 
quantum fluctuations of the valence band. Here we en- 
counter a trilayer generalization of the "orbital" Lamb 
shift, discussed earlier for bilayer graphene. The split- 
ting among the shifted energies {e^} acts as a quantum 
orbital breaking that generally exceeds intrinsic spin or 
valley breaking in scale, and essentially governs the struc- 
ture of the PZM sector. 

The orbital Lamb shift of the PZM Landau levels is 
a "field-theoretic" vacuum effect but is intimately cor- 
related with the Coulomb interaction acting within the 
multiplet. This is because they have to combine to yield 
an electron-hole symmetric spectrum for the PZM multi- 
plet (with only the leading couplings 70 and 71 kept) as a 
whole. In particular, large Coulombic gaps, expected at 
v = and ±3, are essentially given by the energy scale 
~ 2e^ of the orbital Lamb shift. 



The PZM levels get mixed via the Coulomb interaction 
and avoid level crossing, keeping smaller orbital gaps (of 
magnitude ~ 0.5 V c ), as we have seen in Sec. V. Level 
crossing, if present, would enhance the degree of degen- 
eracy and the steps of Hall plateaus would jump ac- 
cordingly. Observations of possible v = ±1,±2,±4, ±5 
quantum Hall states in high-quality samples, such as 
suspended or BN-supported ones, under high magnetic 
fields, if achieved, would be direct evidence for the pres- 
ence of orbital mixing without level crossing. It is also 
possible, in principle, to detect the orbital gaps via cy- 
clotron resonance within the PZM sector. 9,24 

In this paper we have focused on ABC-stacked tri- 
layer graphene. We remark that our analysis and conclu- 
sion cannot simply be carried over to the case of ABA 
trilayers, which lacks a direct link between the K- and 
if'-valley expressions [such as Eq. (2)] and which thus 
requires a separate analysis. 48 
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Appendix A: Derivation of Eq. (21) 

In this appendix we present a proof of the complete- 
ness relation J2 n \9p n \ 2 = e l v I 2 in Eq. (21). A simpler 
version of it is the following: 

00 

f£ n f™ = e^V/a 5 km 

n=0 

for integers k, m > 0, which is verified by use of the 
explicit form of in Eq. (14). We show that Eq. (21) 
is essentially reduced to Eq. (Al). 

Let us first look at Eq. (13) and put the (orthonor- 
mal set of) six eigenvectors of H n for each n e (3, 4, ...) 
into the orthogonal matrix T n — (v ni , v„ 2 , v„ 6 ) with 
I tv I = n and a e (1,2,..., 6). The first row of T n is 

(b$,b$, ...,&1 1 6 ) ) = {bn}}, the second row is {d ( n}}, etc. 
These row vectors also form an orthonormal basis. This 
feature is also true for n G (0, 1, 2), except that T n has a 
smaller rank. 

In taking the product J2 n Sp™ S-'p one mav ^ rs ^ sum 
over TV for each fixed n = \n a \. One thereby en- 
counters inner products of the row vectors such as 

EA$l}{$h = ?> ij and T,Ml){dnl} = 0. The 
remaining sum over \n\ is essentially reduced to for- 
mula (Al) and one eventually finds that J2 n 9p n 9-p = 

e £ 2 p 2 /2 S \k\,\ m \ ( Vfc . which lcadg tQ Eq ^ (21). 
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